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Abstract 

We investigate the assymptotic behaviour of the modeling error in approximate 
deconvolution model in the 3D periodic case, when the order N of deconvolution goes 
Cn| to oo. We consider successively the generalised Helmholz filters of order p and the 

Gaussian filter. For Helmholz filters, we estimate the rate of convergence to zero 
thanks to energy budgets, Gronwall's Lemma and sharp inequalities about Fouriers 
coefficients of the residual stress. We next show why the same analysis does not allow 
to conclude convergence to zero of the error modeling in the case of Gaussian filter, 
leaving open issues. 

0\ MCS Classification : 76D05, 35Q30, 76F65, 76D03 

Key-words : Navier-Stokes equations, Large eddy simulation, Deconvolution models. 

4h 1 Introduction 

Direct Numerical Simulations of flows from the Navier-Stokes Equations (NSE) 

u t + V • (u <g) u) - uAu + Vp = f, 

(1.1) V-u = 0, 
<N u(x,0) = u (x), 

\o 

are accurate only for small Reynold numbers. For large Reynolds numbers, flows are 
turbulent and only means or large scales of velocity and pressure fields might be computed 
thanks to turbulent models. 

Large Eddy Simulation (LES) modeling of turbulent flows aims to apply to the NSE a low 
pass filter specified by a convolution kernel 67, leading to the filtered NSE, written in the 
form 

u t + V- (u®U)-i/Au + Vp = f+ V • S(u,u), 

(1.2) V-u = 0, 

u(x,0) = Uu"(x), 

where u = G * u is the large scale velocity, p = G-kp the large scale pressure, 

(1.3) S(u, u) = u (g> u — u ® u, 

is the subfilter scale stress tensor. A modelisation process aims to seek for suitable ap- 
proximations to S(u, u) in terms of u to close System (1.2), that yields a LES model 
[3, 7, 16]. 
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Figure 1: From Chow et al. 2005 [5]. American Meteorological Society. Reprinted with 
permission. 



Most of LES models are over diffusive and trend to underestimate the energy, creating a 
subfilter scale region (SFS). The total error committed is the sum of the numerical error 
NE and the SFS area [5]. To reduce the SFS area, one uses to apply a deconvolution 
operator to the filter [5, 9, 19, 11, 12]. 

The aim of this paper is to estimate the error modeling in terms of the order of the 
deconvolution denoted by N, in the case of the simplified Bardina's model [1, 10, 4], which 
is based on the approximation 

(1.4) S(u, u) « S(u, u) = u(g>u-n<g>u. 

The approximate deconvolution model (ADM in what follows) is deduced from the sim- 
plified Bardina's model by changing approximation (1.4) in 

(1.5) S(u, u) m §tv(u, u) = u u — Dn(u) <g) Dn(u), 
where the deconvolution operator Dn is such that 

N 

(1.6) D N = -£(I-G) n , 

n=0 

while still noting G the operator associated to the kernel G. We always have §o(u, u) = 
S(u,u), and when ||G|| < l 1 then for a fixed u, 

(1.7) lim Sjv(u,u) =S(u,u). 

Let (un,p n ) be the field calculated from approximation (1.5), that is the solution to the 
system 

d t u N + V • (D N (u N ) D N (\i N )) - vAu N + Vp N = f, 

(1.8) V-ujv = 0, 

uat(0,x) = u^(x), 

if any solution exists. Existence and uniqueness of a solution to System (1.8) was first 
proved in [6] when G is the usual Helmholz filter in the 3D periodic case. More generally, if 
one can prove existence and uniqueness of a solution to system (1.8) for any G that satisfies 

1 the operator norm is based on natural energy spaces the fields belongs to, which will be specified latter 



(1.7), it is expected that the sequence (ujv,Pjv)jVeN converges to (u,p) = (Gu,Gp), for 
some solution (u, p) of the NSE. 

Such convergence results has been proved in [2] in the 3D periodic case, when G = G a:P 
is the generalised Helmholz filter of order p with p > 3/4, where 



after having proved existence and uniqueness of (ujv,Pat)- I n Definition (1.9), Tjj := 
27rZ 3 /L, L > being the size of the computational box, and a > is the filter's width, 
usually of same magnitude of the mesh size in a numerical simulation (see [13] for further 
discussions). 

This yields to consider the error modeling = u — u^r, which goes to zero when N 
goes to infinity. It remains the issue of estimating the rate of convergence in terms of N. 
Staying within the 3D periodic framework and the generalised Helmholz filter of order p 
(p > 3/4), we show in this paper that I? and H 1 norms of £n are of order (p(N+ l)) -1 / 4p , 
(see our main result, Theorem 3.1 below). 

To derive this rate of convergence, we first write the equation satisfied by ejy, by sub- 
stracting (1.8) to (1.2), which yields 

(1.10) d t £N + V • {D N e N ® D N u N ) - vAe N + Vrjv = -V • r jv — V • D N u(g> D N e N , 
where =p — p N , and 

(1.11) tn = u (g) u — -Dyvu <8> D^u 

is the residual stress. By using successively an energy budget procedure and Gronwall's 
Lemma, we get an inequality satisfied by the norms of A 1 I 2 Dn£n where A = G^ 1 (in 
terms of operators), from which we deduce an inequality satisfied by the norms of en 
itself (see Inequality (3.30) below). This inequality highlights the role played by the L 2 
norm of the residual stress. 

The weakness of this method is the regularity assumption that should be imposed on 
the field u, which should be in L^^H 1 ). However, such proceedings are similar to usual 
uniqueness proofs about the NSE, always involving regularity assumptions. 

It remains to estimate the L 2 norm of the residual stress (see Inequality (4.9)). We 
carry out this calculation by using Fourier series expansion and calculations outlined in 
Appendix 7, which, if they use elementary real analysis only, are not straightforward 
and were first speculated thanks to numerical and symbolic computations, before being 
rigorously proved. 

We observe that the rate of convergence slows down as p increases in the range [l,oo[. 
Moreover, the resulting bound goes to a constant that only depends on a and u when p 
goes to infinity and N remains fixed. This is consistent with the idea that more large is 
p, then more smooth are the filtered fields, which should enlarge the SFS area. Therefore, 
one needs high orders of deconvolution to reconstruct well the resolved scale area for large 
values of p. 



2 In terms of operators G a p = (I — a 2p A p ) 1 , where A p denotes the p-Laplacien, and A a p = I — a 2p A p — 
G- 1 



(1.9) 
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Then we consider the popular Gaussian filter, 

(1.12) G«(x) = 

often used in LES. Applying the ADM theory for general abstract filters developed in [17], 
we deduce that the ADM is well-posed in the case of the Gaussian filter. Therefore, one 
may ask if there is convergence of the model to the filtered NSE when N — > oo, and if yes 
what is the convergence rate. 

The theory we develop for Helmholz filters, does not apply to the Gaussian filter, because 
of a too strong convergence of its Fourier modes to zero as the wave number increases, 
although this is not an evidence that the convergence does not hold. 

We argue by approximation in showing that the Gaussian filter can be approximated by 

(2|i |2 \ — m 

when m goes to infinity. We show that our procedure is still valid for this sequence of 
filters, and we derive a bound of order (N + l) _4m fro them. This bound goes to a 
constant depending on a and u when m goes to infinity for a fixed N . Therefore, we 
cannot conclude that the deconvolution process converges to the filtered field (u,p) in 
the case of the Gaussian filter. Because of the strong regularisation effect of this filter, 
we may conjecture that if such a convergence would hold, then it should be very low. 
Therefore, the deconvolution process seems to be not appropriate for the Gaussian filter. 
This remains an open issue. 

The paper is organised as follows. We first fix the mathematical framework and recall the 
results of [2] useful for the continuation of the paper. We next detail how to bound the 
error modeling in terms of the residual stress, whose L 2 norm is then estimated by Fourier 
series expansions. We finally consider the Gaussian Filter by showing how to approximate 
it by the G a 77i's, the error modeling of which being then estimated. The paper finishes by 
a technical appendix including key results to derive estimates about the residual stress. 




2 Mathematical framework 
2.1 Space function 

Throughout the paper, v > and a > are fixed and we stay within the periodic case 
framework. The domain of study is the 3D torus 

(2.1) T 3 = R 3 /T 3 where % := 2vrZ 3 /L, 

for some given L > 0, which is the size of the computational box. All the fields we consider 
have zero mean on T3. Let M s be the vector field space 



(2.2) U s = \ w = ( Wl ,w 2 ,w 3 ) = w k e ik ' x : £ |k| 2s |w k | 2 < 



3 In terms of operators G a , m = ( 1 — ttj— A 
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equipped with the Hermitian structure defined by the inner product and its associated 
norm 



(2.3) (w,v),= E l k ! 2 ^k-vt, ||w|U = E Ik| 2s |w k | 

ker 3 * \ker 3 * 

where 

Vk= (h,k 2 ,k 3 ) eT 3 , \k\ 2 = k\ + kl + kl 
and z* denotes the complex conjugate of z. It can be proved (see [14]) that forall s G 

(2.4) H s is isomorphic to H s (T 3 f, (EI,)' = H_ s , 
and we denote 

(2.5) V(w,v)gH_ s xH s , >,v) s =^w k .v* k 

ker 3 * 

the duality pairing. 

Let H s C EI S be the closed subspace of fields valued in M 3 , characterized by 



Hs = < 



w 



E w k e ikx G M s : Vk G 7%*, w k = w_ k and k • w k = 



On can show (see [14]) that 

(2.6) H s = jw : T 3 M 3 , w G H s (T 3 f, V • w = 0, J wdx = o|, 

2.2 Operators 

2.2.1 Kernel and filter 

The general Helmholz filter w = G Q ^*w is defined by the Fourier Series expansion of the 
kernel G a ^ 

( 2 - 7 ) G «A*) = E G k = j 

ker 3 * ^ 11 

Viewed as an operator, one has G a ^ p = (I— a 2p A 2p ) _1 . Furthermore, a given free divergence 
field w being given, w is solution of the PDE problem 

r v -a 2p A p w + w + Vr = w in T 3 , 

^ ' V-w = inT 3 , 

where the Lagrange multiplier r is constant in this case. 

From now, we write G instead of G a ^ p , and we denote in the same way kernel and operator. 
For all s > 0, G defines an isomorphism, 

M s — > M S+ 2 P 

,2 ' } » a: < w = E w k e lk ' x — > w= E G k w k e tk - X , 
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and we set A = G , characterised by its kernel 

(2.10) A(x) = 4e' k ' x , A k = l + a 2p \k\ 2p . 

ker 3 * 

Notice that if w G H s , then w G H s +2p and the restriction og G to H s , still denoted by 
G is an isomorphism that maps H s onto H S+ 2 P . 

2.2.2 Deconvolution 

Let -Dat denote the deconvolution operator, characterised by the Kernel 

D "= E ( I - G ) n = E^e ik ' x , 

0<n<Af ke% 



where, 



/ 2p 1 1 2p \ '"■ 

% = E I^U =(l + a 2p |k| 2 -W 

(2.11) n=0 V1 ^" |K| ' 



/ a 2p \k\ 2p \ N+1 

PN,p,K-l-\ 1+a2p | k p p 



The following holds [2]: 



(2.12) 1 < Djv.k < N + 1, VkeT 3 , 

l + a 2p |k| 2p 

(2.13) Av, k ) « (7Y + 1) J , for large |k|, 

(2.14) lim D NM = N + l, 

|k|-s>+oo 

(2.15) D NM < (1 + a 2p \k\ 2p ) = A k , Vk G 73, 
where ^4k is defined by (2.10). We deduce from (2.12) and (2.14): 

Lemma 2.1. ^4 reaZ number s > freing given, the operator is a isomorphism over 
M s , such that 1 < ||-Djv|| < N + 1. Morover, the subspace of free divergence field H s is 
stable under the action of D^. □ 

2.3 Former Results 

This section aims to recall results of [2] about the system 



d t u N + V • (D N (u N ) ® D N (u N )) - uAu N + Vp N = f, 

(2.16) V-uat = 0, 

ujv(0,x) = u^(x). 

Throughout the paper, we assume that Uo and f satisfy, 

(2.17) u GH , f GL 2 ([0,T] xT 3 ) 3 , 
and a > is fixed. 
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Definition 2.1 (Regular Weak solution). We say that the couple (u^,p N ) is a "regular 
weak solution" to system (2.16) if and only if the three following items are satisfied: 

1) Regularity 



(2.18) u N G L 2 ([0, T]; H 1+p ) n C([0, T];H p 

(2.19) d t u N G L 2 ([0,T];H ) 

(2.20) p^^OO.T]^ 1 ^)), 

2) Initial data 

(2.21) lim ||ujv(*, •) - uo||h„ = 0, 



3) Weak Formulation 

(2.22) Vvel^O.T]^ 1 ^) 3 ), 
cT r rT 



(2.23) 



/ / d t u N v- / D N (u N ) <g> D N (u N ) : Vv + v \ \ Vujv : Vv 
Jo Jt 3 Jo Jt 3 Jo Jt 3 



3 

+ \ T [ vp N -v= r [ ? v. 

o Jt 3 Jo Jt 3 



□ 



Theorem 2.1. ([2]) Assume p > 3/4. Problem (2.16) has a unique regular weak solution. 
Moreover, when p > 1, 

(2.24) d t u N G L 2 ([0,T],H*-l)> PJV e L 2 ([0,T],iF(T 3 )). 

□ 

Theorem 2.2. There exists a weak dissipative solution to the NSE (1.1) 

(u,p)€ [L 2 ([0,T],Hi)nL 2 ([0,T],H )] x L 5 / 3 ([0, T] x T 3 ) 

such that from the sequence (un,Pn)ngn, one can extract a sub-sequence (still denoted 
(u~N,PN)Nm) such that 

( weakly m L 2 ([0, T], H 1+P (T 3 ) 3 ) D L°°([0, T], H p ), 
u N -> u ^ 

(2.25) [ sJronsfy in L r ([0, T]; H P (T 3 ) 3 ), VI < r < +oo, 



□ 



p w ^p weakly in L 2 ([0, T]; iJ 1 ^) n L 5 / 3 ([0, T]; W /2p,5//3 (T 3 )), 

3 Estimate of the modeling error 
3.1 Regularity assumption and main result 

Let (ujv,Pat) be the solution of Problem (2.16). We assume that the limit (u,p) = (Gu, Gp) 
of (ujv,p~jv)jveN satisfies the regularity assumption 

(3.1) u = 4u£L 4 (Hi). 
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By Sobolev injection Theorem, we deduce 

(3.2) uGL 4 ([0,T] xT 3 ). 
Since (u, p) is solution to the NSE, one has 

(3.3) Ap = -V-(V-(u®u)) + V-f, 
which yields in the periodic case 

(3.4) p€ L 2 {[0,T] x Q), 
and we derive from the NSE, 

(3.5) ^ueI 2 ([0,T],H4 
Our main result is 

Theorem 3.1. Let bn = u — ujv be the error modeling, and assume that (3.1) holds. 
Then we have 

\\e N (t,-)\\ 2 + a 2 P\\e N (t r )\\ 2 p + i; f\\\Vs N (s, + cF\\Ve N (8, -)\\ 2 p )ds < 

Jo 



^ 16Ca 4 L n 



i „^""ii i4(Hl ) 



where C is a universal constant, as a product of Sobolev constants. 4 
3.2 Modeling error and residual stress 

Let £n and tn be the error modeling and the residual stress defined by 

£ N = U — Utv, 

Tjv = U (g) U — DjyU (g) D^U. 

The equation satisfied by ejy is derived by substracting (2.16) to the filtered NSE (1.2). 
Expressing the right hand side in terms of tat, we obtain 



(3.8) d t SN + V • {D N e N (g> D N u N ) - vAsn + Vr^ = -V • t n - V • (D N u <g> D N e N ), 
where r^ =p — p N . 

The aim of this section is to estimate £n in terms of rjy. It adresses A l / 2 D l Jj 2 £is[ rather 
than £jv> since the natural multiplier to get an energy balance from equation (3.8) is 
AD N e N , and formally (d t e N ,AD N e N ) = (d/2dt)||^ 1 /2 jD V 2 £Ar || . Qnce A 1 / 2 D 1 T j 2 s N is 
estimated, we derive bounds for sn (Corollary 3.1 below) by comparing the norms of the 
various operators we consider. I found 4/V 

Theorem 3.2. The following inequality holds: instead of 

\\A^Dfs N{ t, .)\\ 2 + uf Q \\A^s Ni s, OH?* < 

8 e ^H u ll 4 ,4 (Hl ) Cwr^s ||2 ds instead of 

„ e 7 l|TMS ' jll ° dS ' l/i/: this 

/or all N > and t > 0. □ *° be 

checked 



4 For simplicity, we note L 4 (Hi) instead of L 4 ([0, T], Hi) 
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1 I 

Proof. The proof is based on an energy equality satisfied by A2D n en to which one applies 
Gronwall's Lemma. To do so, we use AD^£n as multiplier in the (3.8) satisfied by en 
and we integrate by parts. 

The proof is divided into three steps. In a first one, we check that AD^en is appropriate 
as multiplier to validate the procedure. In a second one, we perform integrations by parts. 
In a last step, we apply usual interpolation inquality to be in order to apply Gronwall's 
Lemma. 

Step 3.i. Consistency of the procedure. We check the regularity of A 1 I 2 D n ' 2 en and 
each factor in equation (3.8) one after each other, beginning with E]y. The regularity 
assumption (3.1) combined with the regularization effect (2.9) of operator G, gives u e 
L 4 ([0, T], Hi + 2p). Therefore, we have at least by (2.19) about u/v's regularity, 

(3.10) e n £ L 2 ([0,T],H 1+p ) C L 2 ([0,T],H 1+p ), 
where Applying Lemma 2.1 combined with (2.9), we get 

(3.11) AD N e N eL 2 ([0,T],Mi- p ). 

We whish to prove now that each factor in equation (3.8) is at least in 

i 2 ([0,T],H P -i) = (L 2 ([0,T],ei_ p ))' 
(see subsection 2.1). To be synthetic, we write things as: 

(2.19) + (3.5) \ f d t E N £ L 2 ([0,T],n ) if 3/4 < p < 1, 



(3 ' 12) + (2.24) d t e N e L 2 ([0,T],H p ^) ifp>l. 

When 3/4 < p < 1, H ^ Ho H p _i, and when p > 1, H p _i ^ H p _i. In all cases, 
(3.13) dt£ N G L 2 ([0,T],H p _i). 

Similarly, 

(2.20) + (3.4) 1 / Vr N a 2 ([0,T]xT 3 ) 3 if3/4<p<l, 
1 } + (2.24) Vr iY GL 2 ([0,T],^- 1 (T 3 ) 3 ) ifp>l, 

that yields 

(3.15) Vr N GL 2 ([0,r],H p _i). 

From the injection Hi )■ Hi, we deduce 



(3.16) (3.10) ^> Ae n £ L 2 ([0, T] 



Furthermore, as (u,p) is a dissipative solution to the NSE, u G L°°([0, T], Ho), therefore 
u € L°°([0,r],H2p), and by lemma 2.1, we get 

(3.17) D N ueL°°{[0,T],W 2p ), 
from which we conclude 

(3.18) < 2 ' 18 > + lem ^!, 7 ]U^ EA . 6 L~([0.rl,H,). 
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Since p > 3/4, we deduce from Sobolev injection Theorem H p > L 4 (Ta) 3 , that yields 
(3.17) + (3.18) 



(3.19) + (2 9) f^V-(D N e N ® D N u N ) E L°° ([0,T}, l ni- 2p _ l , 

Similarly, 



(3.20) V-(D N u®D N e N ) GL oo ([0,T],e 2p _i). 

Finally, u £ L°°([0, T], Ho) combined with (3.2) and properties of G and already 
mentioned, yields 

(3.21) V-r N eL 2 ([0,T],M 2p -i). 
Bringing together all these results, we conclude that when 



A N = d t SN + V • (D N e N <g> D N u N ) - vAe N + Vr N + V • t n + V • (D N u <g> D N e N ) 

then Ajv G L 2 ([0,T],H p _i). Therefore, the duality pairing p _i(Aat, ADn£n)i-p, which 
makes consistent the multiplication of equation (3.8) by AD^bn- In what follows, we 
omit the subscripts when writing duality pairings. 

Step 3.ii. Energy equality. Since all the operators we consider are self adjoint, the 
following holds (see [15]): 

(3.22) (dte N ,AD N e N ) = ^ t \\AlDle N \\l 

(-Ae N ,AD N e N ) = \\A^D^e N \\j. 

Furthermore, since AD^en has zero divergence, (Vr^y AD^sn) = 0. Finally, as the 
operators commute with the differential operators, 

(V • (D N e N D N w N ),AD N e N ) = (A^V ■ {D N s N ® D N w N ),AD N e N ) = 

(3.23) (A- X V • (D N e N ® D N w N ), AD N e N ) = (V • (D N e N ® D N vr N ), D N e N ) = 

((D N w N ■ X7)D N e N ,D N e N ) = 0, 

because -D/v w iV has zero divergence. Finally, arguing as in (3.23) to eliminate the bar in 
the integrals of right hand side, we get 

(3.24) — WA^dIsnWI + v\\A^D\e N \\l = (t n ,VD n s n ) - ((D N e N ■ V)D N u,D N e N ) 



Step 3.iii. Bounds and Gronwal's lemma. We bound each term of the right hand side of 

(3.24) after each other. From Cauchy-Schwarz inequality combined with Young inequality, 
we get 

(3.25) \(r N ,VD N e N )\ < -|M|g + ^H^ejvlli. 
In the same way, by using Ladyzenskaya's inequality [20] we obtain 

\((D N £ N -V)D N u,D N £ N )\<\\D N e N \\ 2 L A\D N u\\i< 

(3.26) i 2 

\\D N £ N \\^\\D N £ N \\l\\D N -a\\i. 
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The symbol of D^G is equal to PN,p,k £ [0, 1] (see (2.11)). Therefore, we have ||Z?Aru||i < 
||u||i. By Young inequality combined with (3.26), we obtain 

(3.27) \((D N e N -V)D N u,D N e N )\ < ^\\u\\f\\D N e N \\ 2 + j\\D N e N \\l 

We deduce from (2.15) that the symbol of Djy is less than the symbol of A l / 2 D l Jj 2 , which 
leads to 

(3.28) \\D N e N \\ <\\A\Dle N \\ , 

regardless of N, Combining (3.24), (3.25), (3.26) and (3.28) yields I found here 

-||r|| 2 + 

(3.29) j\\AlDle N \\l + v\\AlDle N \\l < ^||r||g + h\n\\{\\Al d\s n \\1 V\\J\M\ A Il 

II I I 1 I I 

Inequality (3.9) results from inequality (3.29) thanks to a standard generalisation of Gron- : TO BE 
wall's lemma [8]. □ CHECKED 

Corollary 3.1. The error modeling en satisfies 

\\s N {t,.)\\l + a 2 v\\e N {t,-)\\l + v [\\\Ve N ( S ,-)\\ 2 Q + a 2 P\\Vs N (s,-)\\ 2 p )ds < 
(3-30) Jo 1 t 

/ iiT^( a ,.)n8^ 

v Jo 

for all N > and t>0. □ 
Proof. Let v = ^ v k e ikx G M p . We observe that 



AGAIN 



keT 3 



(3.31) ||^v||g = £ (1 + a 2 P|k| 2 P)|v k | 2 = ||v|| 2 + a 2 ^||v| 



2 

\p- 

ke75 



1/2 

We first take v = sn in (3.31). By using (2.12), which yields the general formal 

1/2 

inequality ||w|| s < \ \D^ w|| s , we deduce the further inequality 

(3.32) \\e N \\ 2 + a 2 v\\e N \\ 2 p < \\A l l 2 D^ 2 e N \\l 

1 /2 

We next take v = diD^ en in (3.31), which yields 

(3.33) \We N \\ 2 + a 2 f||V£Ar|| 2 < WA^D^snW 2 . 
We deduce (3.30) from (3.9) thanks to (3.32) and (3.33). 

□ 



4 Residual stress and rate of convergence 

Now that we have shown that the modeling error ej\r is driven by the L 2 norm of the 
residual stress tn, involving the L 4 (Hi) norm of u, it remains estimate the L 2 norm 
of r^v, which what we aim to carry out in this section. Framework, assumptions and 
notations are those of section 3. 

In what follows, S s denotes the Sobolev constant 5 in the injection H s <^-» L s (T3) 3 . To 
begin with, we show 

5 The constants Si and S1/2 do not depend on L. One can prove that Si < (16 + 3/tt) 1 / 3 , see [14]. 
Unfortunately, we do not know any numerical bound for S1/2, even such a bound may probably be found 
in the litterature 
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Lemma 4.1. The following inequalities holds true 

(4-1) ||T A r||g<2C||u(t,-)||?||u- J D J vu||? /2 , 

(4.2) llu - D N u\\?,~ < - — - — a — /n llull?, 
v. I 1 1 Ml / 2 - (2p(N + 1))V2P M 

where C = SiSu2® D 
Proof. Step 4.i. Proof of (4.1). We write t^v as 

(4.3) TAT = (u-D A ru)Ou + D A ruO(u-D7vu). 

Therefore, combining Holder inequality with 1/3 + 1/6 = 1/2 for conjugation, to the 
Sobolev inequality ||w|| L 6 < Si||w||i, we get 

(4-4) ||t||§ < 25i||u||?||u - D N u\\ 2 L3{j3)3 , 

To estimate ||u — .Divu||2 3 , we use the injection of Wi/ 2 onto L 3 (Ts) 3 to obtain 

(4.5) ||u - D N u\\ 2 L3(j3)3 < S 1/2 \\u - D N u\\ 2 l/2 , 

hence (4.1) by combining (4.4) and (4.5). □ 
Step 4.ii. Proof of (4.2). We deduce from (2.11), 

_ / a 2p| k |2p \2(2V+1) 

(4.6) \\u-D N u\\l /2 =J2{ 1 + JJ k|2 J Ikllukl 2 , 

We apply the technical inequality (7.6) proved in Appendix 7 below, with x = a p |k| p , 
a = 2p(N + 1) > 1, b = 0, which yields 

( ' } U + a 2p |k| 2p y ~ ^(iV + l)' 

We raise both sides of (4.7) to the power 1/p, we multiply the result by |k||uk| 2 and get 
(4-8) , , ' ' |k||u k | 2 <— — — — |k| 2 |u k | 2 , 



l + a 2 P|k| 2 Py 1 11 K| - (2p(N + 1))W 

hence (4.2) from (4.6). 

Corollary 4.1. The following estimate holds 
(4.9) l|r,(V)llg< (2;)(iv 2 ^ 1))1/2p l|u(V)li;. 
for allte [0,7]. 

Inequality (4.9) results from (4.2) combined to (4.1). □ 
Summarizing: (3.30) + (4.9) =>■ 

IMt, 01 IS + » 2p \\e N (t, -)\\ 2 p + i/ j\\\Ve N (s, ■)\\ 2 + a 2 P\\Ve N (s, .)|g)da < 
^ 4 ' 10 ^ 16Ca M ll4 ^||u||4 

11 4 f)^ 3 "l 4 H. 

I/(2 P (iV+l))l/2p l|U ^ 4 (Hl) e 

for all N > and i > 0. 



6 Inequalities (4.1) and (4.2) both hold at any fixed time t £ [0, T], which is not indicated here to reduce 
the notations. 
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5 Case of Gaussian filter 



5.1 Framework 



The Gaussian filter is specified by its kernel, 




(5.2) G(x) = G k e ik - X where G k = 



Let s > and q > s. There exists a constant C be such that 



(5.3) 



Vke7i*, G k |k| 9 <C|k 



Therefore, 



(5.4) 



Vs>0, Vu£i s , Vg>s, GuGH,. 



Let u being given such that Vk S 7^*, |u k | = |k|~ 1-<J ^ (q > 0). Such a vector field u 
belongs to H 9 , but it easy checked that G^ 1 ^ M s for any s. This is why the theory above 
about Helmholz filters fails, since it is based on the fact that G defines an isomorphism 
between H s spaces. 

However, ADM may be considered for the Gaussian filter, and the resulting model yields a 
well posed problem [17]. Moreover, we shall show in what follows that it can be approached 
in some sense, by a sequence of operators which fall within the framework of the theory 
exposed above. 

5.2 Approximation of the Gaussian filter 

we note that for all k £ T% fixed, 



(5.5) 




where 




Let G m denotes the kernel 



(5.6) 




which corresponds to the operator 



still denoted by G, 



(5.7) 




In a sense that needs to be precised, the sequence (G m ) mG ^ converges to G. To be more 
specific, 
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Lemma 5.1. For all k G T3 , 

2 

(5-8) |Gk — G m \.\ < — . 

m 

Proof. We prove in Appendix 7 the technical inequality (7.7) 

Vx>0, Vm>l, fl + — 
\ m 



2 

< — . 



a 2 |k| 2 

We deduce inequality (5.8) in replacing in this inequality x by — — — . □ 

The following corollary is straightforward: 
Corollary 5.1. For all u£i s , 

~ 2 
(5.9) ||Gu - G m u|| s < — ||u|| s . 

m 

In other words, there is weak star convergence of the sequence of operators (G m ) m6 N to 
the Gaussian filter G in H s (s > 0). 

5.3 Powers of the second order filter 

In what follows, we put for m fixed, 



9 a 2 

5.10) n 2 = — — . 

v ; P 24m 

and we denote by H m the m th power of the second order Helmholz operator 

(5.11) H m = (l-^ 2 A)~ m . 

Estimating the error modeling that corresponds to H m yields estimates for the error mod- 
eling that corresponds to G m . The theroy developed above about Helmholz operators 
applies to operator H m . Indeed, let 

(5.12) H mM 



;i + ^ 2 |k 



2\m 



be the symbol of H m . Using the scalar inequality 1 + x m < (1 + x) m < 2 m 1 (l + x m ) for 
positive x, we get 

1 ~ 1 

( 5 - 13 ) nm-1/1 , ..?. m \1,\?.rn.\ - -^"i.k - 



2 m ~ 1 (l + ^ 2m |k| 2m ) ~ m ' ~ 1 + /u 2m |k| 2m ' 

Using results of [2] (section 6), we deduce from (5.13) that the ADM corresponding to 
H m has a unique regular weak solution (un,p n ), in the meaning of Definition 2.1 with 
p = m. Furthermore, this sequence of solution converges to some (u,p) solution of the 
filtered NSE when goes to infinity. Thus we can perform the programme to estimate 
£iv in this case. We next prove. 



7 This estimate is uniform in k, but unfortunately we cannot conclude from this the normal convergence 
of the kernel sequence (G m ) me n because the serie 1/m is not convergent 



14 



Theorem 5.1. Let sn = u — u^r be the error modeling corresponding to H m , and assume 
that (3.1) still holds. Then we have 8 

INivC*, Olio + M 2m ||^(t, Ollm + ^^(HV^^, Olio +M 2r "l|V £ iv( S , Oll^)^ < 
(5 ' 14) Uclm^ 4 

I/(4(iY + l))V2m l|U ^ 4 (Hi) e 



^ll U H 4 ,4 (Hl) 



□ 



Proof. Thanks to (5.13), one can copy line by line proofs of Theorem (3.2) and Corollary 
(3.1) and derive 

\\s N (t,-)\\l + ^ m \\s N (t,-)\\l + u f\\\Vs N (s,-)\\l + ^\\Vs N (s,-)\\ 2 Jds< 

(5-15) J ° o i 4 rt 

-e^ m ^) / ||T^(.,-)llod». 
v Jo 



It remains to estimate ||tjv(s, Olio- Step 4.i in the proof of Lemma 4.1 can be recycled, so 
that (4.5) still holds in this case. Therefore, we only have to bound 



/ ! \2(JV+1) 

(5.16) \\u-D N u\\ 2 1/2 = (l- / 1 + u 2i k i2^ ) |k||u k | 2 , 

ker 3 * v 

where as usual u = X^keT 3 * u ke jkx - We apply the technical inequality (7.2) proved in 
Appendix 7 below, with x = /x|k|, a = 2(N + 1) > 1, m > 1. We obtain 

/ 1 \2(^+l) 

(5 ' 17) V 1 " (1+^lff j ~ (4(iV+l))V2ml k l- 



-. \ 2iV+2 / — 

(5-18) (1- ,, , ) |k||u k | 2 < ^.^"^^ IklVkl 2 , 



We multiply the result by [kHu^l 2 and get 

-y \ 2AM-2 

(I + fl 2 \k\ 2 ) m ) |k| '" KI - (4(iV+l))V2m' 

hence 

2 \[rH[L 
(5.19) u — D/vu L /9 < — — — - u i, 

^ > II iV Ml/2 - ^ N + 1 ^i/ 2m M IU' 

which yields by (4.5), 



(5-20) ||t||S< " " J |u||f, 



|2 ^ 

(4(iV + l))V2r. 



giving (5.14) thanks to (5.15). □ 



3 The constant C below is as in Theorem 3.1, inequality (3.9), see also Lemma 4.1. 
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5.4 Passing to the limit 

From the results of subsection 5.3, we deduce thanks to the relation (5.10) that the ADM 
associated to the filter specified by (5.3), has a unique solution (upf >m ,p N m ) which con- 
verges to some solution (u m ,p m ), of the filtered NSE, by assuming that (u m ,p m ) satisfies 
the regularity assumption (3.1). 

Let £jv,m = u m — UAr m denotes the corresponding error modeling. Thanks to (5.14), we 
obtain 9 

\\e Nim (t r )\\l + a 2m (24m)- m \\e N , m (t,-)\\ 2 m + 



(5.21) 



v I (\\Ve N>m (8,-)\\l + 



o 



a 2m (24m)- m \\Vs N , m (s,-)\\ 2 m )ds < 



70Ca ,. ,, 4 4r||u|| 4 ,, 

11 4 , of 3 " "L 4 (Hi) 



U(4(N+ l))l/2m' 



Without any convergence result about ADM's associated to Gaussian filter (5.1) when iV 
goes to infinity, we cannot consider the corresponding error modeling, and therefore take 
the limit in (5.21) when m goes to infinity. Nervertheless, we observe that for a fixed N, 
the r.h.s of (5.1) converges, as m — > oo, to some C = C(v,a,u,C), which do not depend 
on N. We only can deduce a bound about the sup limit of the terms in the r.h.s. 



6 Conclusions and open problems 

6.1 Typical size of the constants 

The main estimate (3.6) we get in the paper yields the rate of convergence to zero of the 
order modeling in the case of Helmholz filter of order p. The bound involes a constant of 
the form 

(6.1) « = ^iiu|| 4 l4(Hi )^" u " 4l4(Hi) - 

The number of iteration N requiered to reduce substancially the SFS area is driven by the 
size of the constant n. 

This constant involves gradients of the true velocity of the fluid, which may be huge. For 
instance, in some turbulent boundary layer, one may observe flows for which Vu is of 
order 3.10 4 s _1 in layers thick of about 10 _1 m. For such a air layer at 50° (that can be 
considered as incomrpessible) of width and length equal to 1 m, over a time range of 1 s, 
with v w 20.10 -6 m 2 s _1 , we find 

^10 lo28 mV 2 , 

which is a very huge constant. Therefore, even if the resolution would be of order a = 10~ 18 
m, to fully solve such a flow, the number of iteration iV required to substancially reduce 
Sat is so large that the deconvolution algorithm seems not suitable for practical simula- 
tions, which is in contradiction with results of [19], suggesting that very few iterations are 
sufficient to significantly reduce the SFS area. 



5 For the simplicity, we use j^m 1 ^ 2 < 5a instead of (5.10) 
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The rate of convergence as (p(N + l)) _1 / 4 f comes from estimating norms of the residual 
stress ttv involved in the equation for e^, whereas the constant k considered above comes 
from Gronwall's Lemma, which is known to lead to non optimal results. This yields 
the conjecture that the rate of convergence we found is optimal, which is not the case 
of the constant, that might be substantially improved. Furthermore, how the regularity 
assumption u E L 4 (Hi) could be prevented ? 

6.2 Gaussian Filter 

It also remains the issue of convergence of ADM in the case of Gaussian filter. We con- 
jecture that the convergence holds, but in a very weak sense, according to Corollary 5.1, 
a weak sense as yet undefined. 

7 Appendix 

This technical appendix aims at proving a general inequality that has been used in the 
proof of the estimate (4.9). The result is the following. 

Theorem 7.1. The scalar inequality 

i \ a 
1 \ mx 



holds true for any x > 0, a, m > 1. □ 
We consider the LHS function 



h(x) = 1 - t - 



and fixed parameters a,m > 1. 
Its derivative is 



h'(x) = am ( 1 — 



a-1 



1 



x) m J (x + l) m+1 
We apply to h(x) the Lagrange intermediate formula on [0, x] and get 

a-1 

(ip + l) m+1 



(l-- = (x - 0)ti(ip) = xam ( 1 - ' 

V (i + x)"v v ' ^' v (i+vo"v 



for some ip G (0, x). 
The inequality becomes 



i \ a— 1 i 

1 \ 1 mx 



xam 1 — — — —, -r — r-r < 



(1 + ijj) m J (i/j + l) m+1 ~ ^/a 
i.e. (after reducing xm from both sides) 



I \ a_1 1 1 



" ' ' (1 + (l(> + l) m +! ^ 



So now it's enough to prove that 
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; \ a ~ 1 l l 



" ' (l + x) m ) (x + l) m+1 ~ Vfti 



for any x > 0, a, m > 1. 

To easy computations we make the substitution 

1 



6(0,1) 



1 + x 

and the inequality becomes 

a( l_ y m r l ym+ l<^ 

or, after putting a on the RHS 

(1 - y m f- 1 y m+1 < a' 1 - 1 /™ 

for any y £ (0, 1), a, m > 1. 

We denote the LSH above by a(y) = (1 - f)"" 1 y m+1 
Its derivative with respect to y is 

g'(y) = -y m (l - y m ) a ~ 2 ((am + l)y m - (m + 1)) 
We see that the derivative vanishes at 

l/m 



m + 1 



^ am + 1 / 

(\ l/m / \ l/m 

^tl_ J ] th en negative on [ ( j , 1] therefore the max- 



(\ l/m 
) an d is equal to 

/ m + 1 \ l l m / m|l \ a_1 / m+1 \ m / am -m\ a 1 / m + 1 
\am + l/ \ am+lj \am+l/ \ am + 1 / \am + l 

So now we need to show that 

. a-l ■ — 2^ 

am — m 



m + l_\ - <a -l-l/m 



am + 1 / \am + l y 
for a, m > 1. 

Now polish a bit the formula above. In the first term on LHS we simplify m, 



a — 1 \ /m + 1 



a + l/m J \am + l 
In the bottom of the second term on LHS we pull out a, 



< a- 1 " 1 /" 1 



a ~) r'q-i-iMf m + ] ) m <a-^l™ 
a + l/mj \m + l/a/ 



m + 1 
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then cancel a 1 1 / m 



\ a— 1 / , 

a — 1 \ / m + 1 \ m ^ 



a + l/mj \m + l/a 
Then in the second term on LHS we simplify m, 



m+l 



a — I \ a 1 / 1 + 1/m 



a + l/mj \l + l/(ma) 



< 1 



-i \ a— 1 / -i . \ 1+2 

q - 1 \ / 1 + z \ 



Now let z = 1/m G (0, 1] 
The inequality becomes 



v a + z J \1 + z/a, 
for any q > 1, z G (0, 1]. 

We apply the natural log to both sides. We need to show that 

(a - l)/n(a - 1) - (a - l)/n(a + z) + (1 + z)ln{\ + z) - (1 + z)ln(l + z/a) < 
for any a > 1, z G (0, 1]. 
Let 

/(a) = (a - l)/n(a - 1) - (a - l)/n(a + z) + (1 + z)/n(l + z) - (1 + z)ln(l + z/a) 

be the LHS in the inequality above as a function of a. 
The derivative of f (with respect to a) 

f'{a) = ln(a - 1) - ln(z + a) + 

The second derivative is 

„ (z + l)(az-z-a) 

7 { ' a 2 {a-l)(z + a) 

Obviously, since a > 1, z G [0, 1) we have that az — a < 0, so az — z — a < 
therefore 

n, x (z + l)(az - z - a) 
J y ' a 2 (a- l)(z + q) ~ 

We conclude that the first derivative is increasing, therefore 

f'(a) < Urn f'(a) = Urn ( ln(a — 1) — ln(z + a) H — — ) = lira (in ( ] H — — ) = 

a-5>oo a->oc \ a J a^oo \ \Z + a J a J 

Therefore /' is negative, so f is decreasing. It follows that 

f(a) < limf(a) = lira ((a - l)ln(a - 1) - (a - l)ln(a + z) + (1 + z)Zn(l + z) - (1 + z)Zn(l + - j) 
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We know that in general 
therefore, the above limit is zero 



lira xln(x) = 

X->0 



lim(a - l)ln(a - 1) - (a - \)ln(a + z) + (1 + z)ln{\ + z) - (1 + z)ln{\ + -) = 
a— »i a 

= - + (1 + z)Zn(l + z) - (1 + z)Zn(l + 2) = 
We conclude that /(a) < which proves the inequality. □ 
Corollary 7.1. T/ie scalar inequality 



(l + x 2 ) m 7 ~ 2 V2^ 

/ioWs £rae /or any x > 0, a, m > 1. 

In the previous inequality we replace x with x 2 and get 



mx 2 



(7 - 3) I'-OWJ 5 v « 

for any x > 0, a, m > 1. 

Replace in this inequality a with 2a, still keep a > 1 (but works for a > 1/2) 



2a 2 
f ■ > TlTtX 

(7-4) 1- 7, , o W < 



(i + x 2 ) m y - 



for any x > 0, a, m > 1. 

Now extract the square root of both sides 



(l + x 2 ) m y ~ 2 ^/2a 
Remark 7.1. Setting m = 1 in the previous inequality gives 



x 2 \ X . X 



( 7 - 6 ) — o < < 
/or any x > 0, a > 1. 

The following inequality will be used to approximate the Gaussian filter with a power of 
the second order Helmholz filter and calculate the accuracy of this approximation. 

Theorem 7.2. The scalar inequality 

2 



(7.7) \(l + x/n)- n -e- x \< 

is valid for any real x > and any integer n > 1 . 



n 
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It is well-known that as a function of n (and fixed x > 0) the expression 

(l + x/n)- n 

is decreasing and converges to e~ x as n — > oo. (this is elementary calculus, i ommited the 
proof.) 

Therefore, the left hand side in (7.7) can be written as 

|(1 + x/n)- n - e~ X \ = (1 + x/n)~ n - e~ X = e ~^n(l+x/n) _ e -x = e -nln(l+y) _ e ~ny 

where y = x/n > . 

Applying the intermediate value theorem of Lagrange (corresponding to the function £ — > 
e~ n ^) to the last term above we get that 

e -nln(l+y) _ e ~ny = ne ~n^ y _ ^ + y ^ 

for some £ G [ln(l + y),y]. Here we used ln(l + y) < y for y > 0. 
Since e~ n ^ < e - nln ( l +y) we further have that 



(7.8) e-^+v) - e -ny < ne -nln(l+y) ^ _ ^ + y y = n(1 + y yn^ y _ ^ + y)) 

for any real y > and integer n > 1 

The term y — Zn(l + y) appearing in the last term in the inequality above is estimated as 

0<y-ln(l + y) < y - 

for any real y > 

Going back to inequality (7.8) we finaly have 

,.2 

e -nln{X+y) _ e ~ny < ^ + ^-nfL 

We replace y = x/n and obtain 



l + -)""-e-<n(l + ^- n -^=x 2 (l + ^- n ^- 

But, as pointed out before, for any fixed x the function n — >• (1 + x/n) _n is decreasing, so 
we have that for n > 2 

x\~ n / x\~ 2 4 
1 + -) <(1 + ^1 < 



nJ V 2/ 1 + x 2 

Therefore, for n > 2 

x\~ n 2 x\~ n 1 „ 4x 2 12 

1 + - -e~ x <x 2 1 + - — < s — <- 

nJ V n/ 2n 1 + x^ 2ra n 

For n = 1 the left hand side of (7.7) becomes 

(1 + x)- 1 - e~ x < (1 + x)" 1 + e" x < 2 
for any x > 0, so the inequality (7.7) is valid for n = 1 too. 
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